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Abstract 

We consider three-loop radiative corrections to the Lamb shift and hyperfine splitting. Correc- 
tions of order a 3 (Za) 5 m are the largest still unknown contributions to the Lamb shift in hydrogen. 
We calculated radiative corrections to the Lamb shift and hyperfine splitting generated by the 
diagrams with insertions of one radiative photon and electron polarization loops in the graphs 
with two external photons. We also obtained corrections generated by the gauge invariant sets 
of diagrams with two reducible radiative photon insertions in the electron line and polarization 
operator insertion in one of the radiative photons, and diagrams with two reducible radiative pho- 
ton insertions in the electron line and polarization operator insertion in one the external photons. 
Corrections to the Lamb shift and hyperfine splitting generated by the diagrams with insertions 
of the three-loop one-particle reducible diagrams with radiative photons in the electron line are 
calculated in the Yennie gauge. 
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I. INTRODUCTION 



The largest still unknown contributions to the Lamb shift in hydrogen are corrections of 
order a 3 (Za) 5 m. The magnitude of these corrections can be easily estimated multiplying 
the corrections of order a 2 (Za) 5 m (see, e.g., [1,2]) by an extra factor a/ir. The contributions 
of order a 3 (Za) 5 m turn out be as large as 1 kHz for the 15* state in hydrogen. Below we 
present our recent results for three-loop radiative corrections to the Lamb shift and hyperfme 
splitting. 

Nonrecoil corrections of order a 3 (Za) 5 m to the Lamb shift and corrections of order 
a 3 (Za)Ep to hyperfme splitting are generated by three-loop radiative insertions in the 
skeleton diagram in Fig. [1] Respective corrections of lower orders in a generated by one- 
and two- loop radiative insertions are already well known (see, e.g., review 0]). The crucial 
observation, which greatly facilitates further calculations, is that the scattering approxima- 
tion is adequate for calculation of all corrections of order a n (Za) 5 m and a n {Za)Ep (see, 
e.g., a detailed proof in Q). One may easily understand the physical reasons which lead to 
this conclusion. Consider the matrix elements of the skeleton diagram in Fig. [T] with the on 
shell external electron lines calculated between the free electron spinors, and multiplied by 
the square of the Schrodinger-Coulomb wave function at the origin. They are described by 
the infrared divergent integral 

16(Za) 5 /m r \3 f°° dk r 

m — S l0 , (1) 



7m 3 V m ) J fc 4 
in the case of the Lamb shift, and by the infrared divergent integral 1 



in the case of hyperfme splitting. In these integrals k is the dimensionless momentum of the 
exchanged photons measured in the units of the electron mass. 



FIG. 1: Skeleton two-photon diagram 

Let us consider radiative insertions in the skeleton two-photon diagram in Fig. [Q Account of 
these corrections effectively leads to insertion of an additional factor L(k) in the divergent 
integrals above, and while this factor has at most a logarithmic asymptotic behavior at 



We define the Fermi energy Ep as 

E F = l ±zW^-(l + a,) (^YchRoo, (2) 
3 M V to / 

where m is the electron mass, M is the muon mass, m r is the reduced mass, a is the fine structure 
constant, c is the velocity of light, h is the Planck constant, R x is the Rydberg constant, a M is the muon 
anomalous magnetic moment, and Z is the nucleus charge in terms of the electron charge [Z = 1 for 
hydrogen and muonium). 



2 



large momenta and does not spoil the ultraviolet convergence of the integrals, in the low- 
momentum region it behaves as L(k) ~ k 2 (again up to logarithmic factors), and improves 
the low-frequency behavior of the integrand. However, the integral for the Lamb shift 
is sometimes still divergent after inclusion of the radiative corrections because the two- 
photon-exchange diagram, even with radiative corrections, contains a contribution of the 
previous order in Za. This spurious contribution should be removed by subtracting the 
leading low-momentum term from L(k)/k A . The result of such subtraction is a convergent 
integral, where the low integration momenta (of atomic order mZa) in the exchange loops 
are suppressed, and the effective loop integration momenta are of order m. Then it is clear 
that small virtuality of the external electron lines would lead to an additional suppression 
of the matrix element under consideration, and it is sufficient to consider the diagrams only 
with on-mass-shell external momenta for calculation of the contributions to the energy shifts. 
As an additional bonus of this approach one does not need to worry about the ultraviolet 
divergence of the one-loop radiative corrections. The subtraction automatically eliminates 
any ultraviolet divergent terms and the result is both ultraviolet and infrared finite. 



II. DIAGRAMS WITH THREE ONE-LOOP ELECTRON VACUUM POLARIZA- 
TIONS 

A. Lamb Shift 

Each polarization loop in the diagrams in Fig. [2] corresponds to insertion of the vacuum 
polarization operator (a/n)k 2 Ii e jEf in the Lamb shift skeleton integral in eq.([T]), and the 
contribution to the Lamb shift generated by the diagrams in Fig. [2] has the form [5J 



(^Y m f dkk 2 l\ e = - 0. 021 458 (1) Q (^Y m. (4) 

^ m ' Jo 11 n \m / 



Numerically for the IS" level in hydrogen we obtain 5E^ = -0.002 16 kHz. 



FIG. 2: Three one-loop polarizations 



B. Hyperfine Splitting 

We obtain the expression for the radiative correction to hyperfine splitting generated by 
the diagrams in Fig. [2] inserting the polarization loops in the skeleton integral in eq. ([3]) ^ 



5E 



= f = 2- 3 (4)2W E F . (5) 

' ' J 



Numerically for the ground state in muonium the correction is 6E$ FS = 0. 003 29 kHz. 
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III. DIAGRAMS WITH TWO-LOOP AND ONE-LOOP ELECTRON VACUUM 
POLARIZATIONS 



A. Lamb Shift 

The integral for the diagrams in Fig. [3] is obtained from the skeleton integral in eq. (JTJ) 
by insertion of the one-loop vacuum polarization (a/ir)k 2 Ii e , and the two-loop vacuum 
polarization (a/n) 2 k 2 l2e |5j 

5E { ? = M m \ dfefie/ae = - 0.390 152 (7) \ ' ( — ) m. 6 

(2) 

Numerically for the IS" level in hydrogen we obtain 5E K L > = -0.039 21 kHz. 



FIG. 3: One- and two-loop polarizations 



B. Hyperfine Splitting 

The expression for the hyperfine splitting generated by the diagrams in Fig. [3] has the 
form [B| 

SE f F3 = «S2f2> B, r dkk%J 2e = 3. 559 9 (2) ^-E F . (7) 

71 Tl Jo 71 

Numerically for the ground state in muonium we obtain 5E£ J FS = 0.004 56 kHz. 

IV. DIAGRAMS WITH THREE-LOOP ELECTRON VACUUM POLARIZATION 
A. Lamb Shift 

For calculation of the correction generated by the diagrams in Fig. H] we need the three- 
loop vacuum polarization operator (a/7r) 3 k 2 I^ e . Seven leading terms both in the low- and 
high-momentum asymptotic expansions of over the powers of the momentum were calcu- 
lated analytically in p, 0, H, @, We adjusted these results for the case of the momentum 



renormalization scheme used in QED, and constructed an interpolation which approximates 
the three-loop polarization operator for all Euclidean momenta. 

The skeleton integral in eq. (pQ) remains infrared divergent even after insertion of the three- 
loop vacuum polarization since /3 e (0) ^ 0. This linear infrared divergence is effectively cut 
off at the characteristic atomic scale mZa if we restore finite virtualities of the external 
electron lines. As mentioned in the Introduction, such infrared divergence lowers the power 
of the factor Za, and respective would be divergent contribution turns out to be of order 
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a 3 (Za) 4 . We carry out the subtraction of the leading low-frequency asymptote of the polar- 
ization operator insertion, which corresponds to the subtraction of the leading low-frequency 
asymtote in the integrand for the contribution to the energy shift h e {k) = he{k) — l3 e (0), 
and insert the subtracted expression in the formula for the Lamb shift in eq.([T]). Then the 
contribution to the energy shift has the form j^] 

c , 3) 2>2a 3 (Zaf /m r \3 f°° dk ~ . . a 3 (Za) 5 /m r \3 

*4 3) = -z-r^ (-? ) m / ^h e = 1.015 88 (5) -L-J- -f m. (8) 



7r 4 n 3 V m J J Q k 2 n 2 n 3 \ m 

Numerically for the 15* level in hydrogen we obtain 5EY> = 0.102 10 kHz. 



FIG. 4: Three-loop polarizations 



B. Hyperfine Splitting 

There is no problem of infrared divergence for the radiative correction to hyperfine split- 
ting generated by the three-loop polarization insertions in Fig. H] (^] 

= 1^1 Ef I" dk h . = 1.647 9 (5) ^P-E F . (9) 
ft n Jo 71 

(3) 

Numerically for the ground state in muonium we obtain 5Efi J FS = 0.002 11 kHz. 

V. DIAGRAMS WITH ONE-LOOP ELECTRON FACTOR AND TWO ONE- 
LOOP ELECTRON VACUUM POLARIZATIONS 

A. Lamb Shift 

To calculate the correction of order a 3 (Za) 5 generated by the gauge invariant set of 
diagrams in Fig. [5] we need a new element, namely, the gauge invariant electron factor L^ik) 
in Fig. [6] which describes all possible insertions of the radiative photon in the electron line 
with two external pho tons. An explicit expression for this electron factor was obtained in 



different forms in [111 . Il2l . Il3l . |l4|. Inserting in the skeleton integral in eq.([T]) the electron 
factor (a/ir)k 2 LL(k), one-loop polarization operator squared and the multiplicity factor 3 
we obtain the radiative correction in the form 

SE m^j^W ( ^y m f°° dkk , Ldk)Il (10) 



o 



It is easy to check explicitly that this integral is both ultraviolet and infrared finite. After 
numerical calculations we obtain |5l 
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6E^= 0.0773 (4) (^Y m. 

For the 15* level in hydrogen this contribution is SE^ = 0.007 77 (4) kHz. 




FIG. 5: One-loop electron factor and two one-loop polarizations 




FIG. 6: One-loop electron factor 



B. Hyperfine Splitting 

We calculate the contribution to hyperfine splitting generated by the diagrams in Fig. 
[5] using an explicit expression for the electron factor Lhfs like in the case of the Lamb 
shift above. This electron factor was obtained in [151 ] . Inserting the electron factor 
(a/ir)k 2 LHFs(k) in the skeleton integral in eq.Q we obtain the radiative correction in 
the form [5j 

5E f FS = 2 -^^E F dkk'L HFS (k)ll = - 3.487 2 (2) ^±E F . (12) 
Numerically for the ground state in muonium this contribution is 6E$ FS = -0.004 47 kHz. 

VI. DIAGRAMS WITH ONE-LOOP ELECTRON FACTOR AND TWO-LOOP 
ELECTRON VACUUM POLARIZATION 

A. Lamb Shift 

An integral representation for the correction generated by the diagrams in Fig. [7] is 
obtained from the skeleton integral in eq. ([1]) in the standard way [E| 



l(5) 32a 3 {Za) 5 (m r \z _ f°° JUT n ^ T , ni ^ a 3 (Zaf fm r \\ 



SE ? = ~ S^UfJ m l *W,.= 2.1913(4)^^) m . (13) 



Numerically for the IS* level in hydrogen this contribution is 8E^ = 0.220 24 (4) kHz 
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FIG. 7: One-loop electron factor and two-loop polarization 



B. Hyperfine Splitting 

The radiative correction to hyperfine splitting generated by the diagrams in Fig. [7| has 
the form \^ 

SE W = \Qa\Za) Ep r°° dkk 2 LHFs{k)he = _ 46g0 9 (1) ^1M Ef . (14) 

7T 71 Jo IT 

(5) 

Numerically for the ground state in muonium we obtain 6E% J FS = -0.006 00 kHz. 



VII. DIAGRAMS WITH ONE-LOOP POLARIZATION INSERTIONS IN THE 
ELECTRON FACTOR AND IN THE EXTERN ALAINED PHOTON 

A. Lamb Shift 

The contribution to the Lamb shift generated by the diagrams in Fig. [8] is similar to 
the contribution in eq. fll3l) . the only difference is that now we consider a radiatively cor- 
rected electron factor in Fig. [9] and a one-loop polarization insertion in the external photon. 
Insertions in the skeleton integral in eq.([T]) lead to the expression j^] 

*#--^(sy~r**»w»- 0.03 7 3 6 (1 )^eV ( 15 ) 

where the the electron factor Z^ 2 ' 1 ) with one-loop polarization insertion in Fig. was obtained 



m 



13| . Numerically for the 15* level in hydrogen this contribution is SE^ = 0.003 75 kHz. 



FIG. 8: One-loop polarization insertions in the electron factor and external photon 



rOn /X, /X, 

i i = i r + 2 J —r L — t + ~i — L - 



FIG. 9: One-loop polarization insertions in the electron factor 
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B. Hyperfine Splitting 

We obtain the radiative correction to hyperfine splitting generated by the diagrams in 
Fig. [H] inserting the radiatively corrected electron factor (a/Tr)k 2 L^'p S (k) [l6| in Fig. in 
the skeleton integral in eq.® [B| 

SEP FS = l -^^E F ^ dkk>L™(k)I u = -0.533 3 (5) ^^E F . (16) 

7T Tl Jq IT 

Numerically for the ground state in muonium this contribution is 6E$ FS = -0.000 68 kHz. 

VIII. DIAGRAMS WITH TWO ONE-LOOP POLARIZATION INSERTIONS IN 
THE ELECTRON FACTOR 

A. Lamb Shift 

The contribution to the Lamb shift generated by the diagrams in Fig. [10] is similar to the 
correction generated by the one-loop polarization insertion in the electron factor calculated 
in [ljjj]. The explicit expression for this correction contains the electron factor with two 
one-loop polarization insertions L^ 3,1 ^(fc) in Fig. [11] j^] 

SEl = M-{m) m l dk P = -°- 012610(3) ^^tJ 

(17) 

(7) 

Numerically for the IS" level in hydrogen we obtain 5E}"' = -0.001 27 kHz. 



FIG. 10: One-loop polarization insertions in the electron factor 



i i = i r + 2 - C -7 J — r + r i 

FIG. 11: Two one-loop polarization insertions in the electron factor 



B. Hyperfine Splitting 

The contribution to hyperfine splitting generated by the diagrams in Fig. [10] is similar to 
the correction generated by the one-loop polarization insertion in the electron factor which 
was calculated in 1^|. The explicit expression for this correction has the form j^] 



SE g = So^o) Ep r°° dkL (^ {k) = _ Q 309 Q5 (7) ^1 Ef _ (18) 
n n Jo ' 71 
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(7) 

Numerically for the ground state in muonium we obtain 5E% FS = -0.000 40 kHz. 

IX. DIAGRAMS WITH TWO-LOOP POLARIZATION INSERTION IN THE 
ELECTRON FACTOR 

A. Lamb Shift 

The contribution to the Lamb shift generated by the diagrams in Fig. [12] is similar to the 
correction generated by the one-loop polarization insertion in the electron factor which was 



calculated in [13]. The explicit expression for this correction 



differs from the respective expression in 13[ only due to the difference between the electron 
factor with one- loop polarization insertion L L ' (k) in Fig. [9] and the electron factor with the 
two-loop polarization insertion L L ' (k) in Fig. [13j After numerical calculations we obtain 
1 

5E$> = -0.245 71 (7) a3( ff (^X m. (20) 



7r 2 n 3 V m 

(8) 

For the 15 level in hydrogen the correction is 5E L = —0.024 70 kHz 



FIG. 12: Two-loop polarization insertions in the electron factor 

+ 2 r ? 1 + 4 — r + 2 — V 
FIG. 13: Two-loop polarization insertions in the electron factor 



B. Hyperfine Splitting 

The contribution to hyperfine splitting generated by the diagrams in Fig. [12] is similar to 
the correction generated by the one-loop polarization insertion in the electron factor which 
was calculated in T^. The explicit expression for this correction has the form j^] 

6EP FS = ^ dkL%%(k) = -0.123 9 (6) (21) 
Numerically for the ground state in muonium we obtain 8E%' FS = -0.000 16 kHz. 
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X. DIAGRAMS WITH FACTORIZED RADIATIVE PHOTONS AND POLAR- 
IZATION OPERATOR INSERTIONS IN ONE THE RADIATIVE PHOTONS 



A. Lamb shift 

We use representations 



A M (fc, t) = ^ { A T k 2 ln + B rltM {p -k-m) + C TPfl {p -k-m) + E T a^k x } , (22) 

£(p — k,r) = £-(p -k-mf \M Tl + (p - k)M r2 \ , 

for one-loop vertex and mass operator with massive photon, and similar representations 
A At (/c) and E(p — fc) for one-loop vertex and mass operator with massless photon to obtain 
an explicit expression for the contribution to the Lamb shift generated by the diagrams in 
Fig.[H 

5E W = Aa?{Zaf /"M 3 m f 1 dv V i 1 s) /• rfx f°° dk[{-A+B+C-E+M x ){Ar-M^) 
71 n \ m / Jq 1 — v J Jq 

+ (-A T + B T + C T -E T + M Tl )(A - M 2 )} = -0.139 97 (1) ^f^' ( — Y m. 
Numerically for the IS" level in hydrogen this correction is SE^ = —0.014 07 kHz. 



+ 2 > > + 2 > > + 2 



FIG. 14: One-loop polarization insertions in one of the radiative photons 
B. Hyperfine Splitting 

The contribution to hyperfine splitting generated by the diagrams in Fig. [TH has the form 

SE p = 2a\za) Ef r 1 dv v (* 3) r dx r°° dk ^ 2]g 2 AAr _ 2 ^ AEr + ArE) (23) 

^n 6 Jq 1 - v 2 J J 

-2(BB T -BE T -EB T + EE T )+C(E T -B T ) + C T (E-B)-2k 2 AM T2 + (2B + C)(-M Tl + M T2 ) 

+2EM Tl -2k 2 A T M 2 +(2B T +C T )(-M 1 +M 2 )+2E T M 1 -2(M 1 -M 2 )(M Tl -M r2 )+2k 2 M 2 M T ^ 

, , a 3 (Za) 
= 0.051 44 4 ; Ep. 

Numerically for the ground state in muonium we obtain 5E$ FS = 0.000 07 kHz. 
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XI. DIAGRAMS WITH FACTORIZED RADIATIVE PHOTONS AND POLAR- 
IZATION OPERATOR INSERTIONS IN ONE THE EXTERNAL PHOTONS 

A. Lamb shift 

The contribution to the Lamb shift generated by the diagrams in Fig. [15] in terms of 
one-loop radiative corrections has the form 

+(-A + B + C-E + M 1 )(A- M 2 )\ = -0. 006 25 (1) \ J m. 
Numerically for the 15* level in hydrogen we obtain 5e£ 0) = -0.000 63 kHz. 

+ 2 O 5 +2 O +2 



FIG. 15: One-loop polarization insertions in one of the external photons 



B. Hyperfine Splitting 

The contribution to hyperfine splitting generated by the diagrams in Fig. [15] has the form 

= jf * r "SfeS / «^ - ,AE (24, 

-2{B 2 - 2BE + E 2 ) + 2C(E - B) - 2k 2 AM 2 + (2B + C){-M l + M 2 ) 
+2EM 1 - 2k 2 AM 2 + (2B + C)(-M x + M 2 ) + 2EM 1 - 2(M l - M 2 f + 2k 2 M 2 2 ] 

= 0.024 8 (5) ^E F . 

Numerically for the ground state in muonium we obtain 5E%% = 0.000 03 kHz. 

XII. DISCUSSION OF RESULTS 

We calculated three-loop radiative corrections to the Lamb shift of order a 3 (Za) 5 m and 
three-loop radiative corrections to hyperfine splitting of order a 3 (Za)Ep generated by the 
gauge invariant sets of diagrams in Figs. El El HI El El [HE CGI EU and [HI Collecting all 
contributions above we obtain 

i=10 3/ „ 

5E L = Y 8Ef = 2.5057 - \ 7 m. (25) 
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Numerically for the IS* level in hydrogen this contribution is 8El = 0.252 kHz. 
For hyperfme splitting we obtain 



5E HFS = V 6E% S = -1.282 ^$E F . (26) 

8=1 

Numerically for the ground state in muonium the contribution is 5Ehfs — —0.001 7 kHz. 

S-> Si Sp £2 S~t 

} I + 2 ? ? + ? ? 



FIG. 16: Reducible three-loop diagrams 

We also calculated corrections of order ot'iZafm to the Lamb shift and corrections of 
order a 3 (Za)E F to hyperfme splitting generated by the diagrams in Fig. [161 These factorized 
contributions are not gauge invariant and calculations we made in the Yennie gauge. We 
obtained [13] 

AE L = -5. 321 93 (1) m, (27) 

or AEl = —0. 535 kHz for the IS" level in hydrogen. 



Respective contribution to hyperfme splitting is 17 



AEhfs = 0. 104 23 (1) a3( f Q) E F , (28) 

or AEhfs — 0. 000 13 kHz for the ground state in muonium. 

Work on calculation of the remaining three-loop contributions to the Lamb shift and 
hyperfme splitting is now in progress. 
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